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Abstract 

Tests for break points detection in the law of random vectors have been proposed in several papers. Nevertheless, they 
have often little power for alternatives involving a change in the dependence between components of vectors. Specific tests 
for detection of a change in the copula of random vectors have also been proposed in recent papers, but they do not allow 
to conclude of a change in the dependence structure without condition that the margins are constant. The goal of this 
article is to propose a test for detection of a break in the copula when changes in marginal distribution occurs at known 
instants. The performances of this test are illustrated by Monte Carlo simulations. 
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1. Introduction 


Let A be a d-dimensional random vector (d > 2), with cumulative distribution function (c.d.f.) F and marginal 
cum ulative distr ibution functions (m.c.d.f.s) Fi,... ,Fd. When the m.c.d.f.s Fi,... ,Fd are continuous, Sklar’s Theorem 
(see Sklar . 1959ll allows us to say that there exists a unique function C called copula, characterizing the dependence of 
random vector X, such that F can be written as: 

F{x) = CiFiixi),...,Fd{Xd)), xgR‘^. ( 1 ) 

Let Xi,... ,Xn be d-dimensional observations. The purpose of change-points detection is to test the hypothesis 

Flo ■ 3F such as Xi ,..., A„ have c.d.f. F, (2) 

against -■ido- The equation (IT|) involves that Hq can be rewritten as Hq = idp.m H ido.o with 


ido,m : 3 Fi, .. .Fd such as Xi ,..., A„ have m.c.d.f. Fi,... Fd, 
Hq^c ■ dC such as Ai,..., A„ have copula C. 


( 3 ) 

( 4 ) 


A change either in the copula of random vectors or in one of m.c.d.f.s implies the rejection of the null hy pothesis Hp . 
Many non-parametric tests for Hp based on empirical processes are present in the literature; see for example iBail ()l994ll : 
Csorgd and Horvhdil ( 1997 ): Inoud ( 200ll) . These tests a re not very sensitive to detect a change in the copula which leaves 
the m.c.d.f.s unchanged. This conclusion is highlighted in iHolmes et al.l (|2013L Section 2) through Monte Carlo simulations. 


Non-parametric tests for break detection, sensitiv e to changes in the c opula of observations and based on the two-sided 
sequential empirical copula process are considered in iBiicher et al.l (120141 ). These tests do not allow to conclude in favour 
of ~'Ho c without condition on the constancy of m.c.d.f.s. In many situations, see for example section 5.2, a specific event 
can lead to changes in the marginal cumulative distributions. The question then becomes whether the specific event 
changes the copula or not. The aim of this paper is to propose a test to detect a change in the dependence structure of 
random vectors, sensitive to changes in copula of observations and adapted in the case of alternative hypotheses involving 
abrupt changes in m.c.d.f.s. 


This paper is organized as follows. The procedure to test the null hypothesis of a break in c.d.f. when a change in 
the m.c.d.f.s occurs is presented in Section 2. An adaptation of results of Section 2 when multiple changes in m.c.d.f.s 
occur is described in Section 3. Section 4 contains the results of Monte Carlo simulations. Finally, Section 5 reports brief 
discussions about the case of a/pha-mixing observations and presents an illustration on a specific situation. 
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2. Break detection in the copula when a break time in the m.c.d.f.s is known 


In the sequel, the weak convergence, denoted by must be understood as being the weak convergence in the sense 
of Definition 1.3.3 in van der Vaart and Wellner ( 20001) . For a set T, £°°{T) denotes the space of bounded real-valued 
functions on T equipped with the uniform metric. 

Let JCi, ..., Xn be d-dimensional rand om vectors ( d > 2 ) and consider for 1 < k < I < n the empirical copula Ck-.i of 
the sub-sample Xk,.. .,Xi as suggested in IPeheuvelsI ( 19791 ): 


l d 


Ckdi'u) — 7 _ 7 I 1 ^ %)) ^ [ 0 ; 1 ] 


I — k + \ 

i=kj=l 

where for j = 1,..., d, Fk-.i,j is the empirical cumulative distribution function (e.c.d.f.) of sample Xkj^ ..., Xij: 

1 ^ 

Fk:l,j{x) = ^ ^ 'y ^ — x), X G R. 


( 5 ) 


i—k 


In lBiicher et al.l (|2014D . the following Cramer-von Mises’s type statistic to test -'Flo is suggested: 

Sn= sup \/nA„(s, 1)A„(0, s) / {CiqnsJ («) - C'LnsJ-Hl:n(M)}^Ci,„('u), 
sefo.n dro.iF 


sG[0,l] 

where Xn{s, t) = {[nt\ — [nsj )/n, s < t € [0,1]. 


( 6 ) 


( 7 ) 


Monte Carlo simulations (see section 5 of lBiicher et ah . 20141) hi ghlighted that a strategy of boostraping with inde- 
pendent or dependent multipliers according to the observations (see iBiicher and Koiadinovid 2015 : Bucher et al.l . [20l3) 
of the statistic Sn leads to very good performances in term of powers for alternatives hypotheses that involve a change in 
copula which leave the m.c.d.f.s unchanged. 

Let us suppose that it exists a break time m = [nb\ in m.c.d.f.s, b G (0,1) known. We propose a test for = 
i7o,c n Fli^rn, where Hq c is defined in (l4|) and Fil m is defined by: 


Hi m ■ 37^1 ,.. .Fd and F^, ■ ■ ■ F^ such that 


Xi,..., Xm have m.c.d.f. Fi,... 
Xm+i, ■■ ■,Xn have m.c.d.f. F{,... F^. 


( 8 ) 


Note that we do not suppose that F^,. ■ -F^ are necessarily different from Fi,... F^. In other words, we do not assume a 
change in the m.c.d.f.s. However we suppose that if there is a change in the m.c.d.f.s, it is a unique and abrupt change at 
time TO. 

Let Xi,. ■ ■, Xn be d-dimensional random vectors with unknown copula C, such that Xi^ ■ ■, Xm have m.c.d.f.s 
Fi,. ■ ■ ,Fd and Xm+i, ■ • ■, Xn have m.c.d.f.s F(,..., F^, where Fi,...,F^, F(,...,F^ are unknown and the break time 
TO = [nb \, & G (0,1) is known. 


For i G {1,..., n}, let us consider the random vectors Ui^m defined by 

Fli ra — 


(Fi(X,iFd(X,d)) 

(F((X,i),..., F’diX^d)) iG{m + l,...,n}. 


(9) 


Note that the vectors Ui^m, f = 1,..., n have C for c.d.f. For i G {1,..., n}, let Ulm defined by: 

jyi-n f {Fl:rn,l{Xii'j ^ ; Fi-nri,diXid')') i G {1, . . . , to} 

\iFm+l-.nA^ii)y ■ ■ yPm+l-.n,d{Xid)) i G {m + 1, . ■ ■ ,n}, 

where for all 1 < A: < ^ < n and j = 1,..., d Fk-i j is the empirical c.d.f. of Xkj ,..., Xij as defined in ([6]). The vectors 
, n can be seen as pseudo-observations of copula C. An estimator of C is given by the empirical distribution 
of rri-” r/i'" • 

n 

Ci:„,„7(«) = -Vl(f//;((<«), «G[0,1]‘'. 

n 

i—l 
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This estimator can be rewritten for u G [0,1]^ by: 


^ / \ n — m 

n n 


where for any subsample , Xi, 1 < k < I < n, Ckd is the empirical c.d.f. of random vectors , Uj^'^ defined 

in m, and for i = k,... ,l 

ay = iFk■.lAx^l),.■.,Fk,lA^^l)). 

More particularly, for 1 < fc < ^ < n, m S {1,..., n — 1} and u G [0, l]'^, 


(7 ('mI — I l-k+l ^k-.miu) + i_f.j^iCm+ld{u) m G [k,l], 
Ck-.iMA \ Caiu) m^[kj]. 


For a subsample Xk ..., Xi, 1 < k < I < n consider the following pseudo-observations of copula C: 

{j {Ff^-Yn^l ); ■ • ■ 5 )) \U I] 

\{Fm+l■.l,l{XA,■■■,Fm+l■.lAX^l)) iG{m+l,...,l} 

{Fk:i^i{XA, ■ ■ ■, Fk-.iA^A) [k,l]i 

Fkd defined in ([6]). Then Ckd,m defined in (fT())l is the empirical cumulative distribution of ,..., UiA- 
The corresponding two-sided sequential empirical copula process is defined by 


( 10 ) 


C„,m(s,t,M) = Vn\Asd){Cyns\+i:lnt\,AA “ C'(«)}: {s,t,u) G A X [0, l]^', 

[nt] 

V ^ - I I , T 

I— [nsj +1 


A = {(s,t) G [0, l]^|s < t}. The test statistic proposed in this paper is based on the process ^n,m, defined by 

lll^n,m('55 ^) — (0, s) An('S, ("u) f-^[nsj(‘^) } ? [0:^] 

Note that V)n,m can be rewritten as 

Dn,m(s, u) = A„(s, 1)C„,,„(0, S, u) - A„(0, s)Cn,As^ 1- “)> (s, «) S [0, 1]^^+^ 

Similarly to S'„ defined in we consider the Cramer-von Mises statistic 


Sn,m — sup / '[ID)yj^m(s, "u)} ^l:n,m(A 
sG[0,l] J[0,1]‘^ 

1 


= .max -J2{^n,m(k/n,UlA)A. 
l<k<n—l n ^' 


( 11 ) 


( 12 ) 


The asymptotic behaviour of the empirical process Dn,m is given o n Proposition proved in Appendix A. The result 
is obtained under the following non-restrictive condition, proposed in ISeger i (l2012ll : 


Condition 2.1. For any j G the partial derivatives Cj = dC/duj exist and are continuous on = {it G 

[0,l]^u, G (0,1)}. 

Proposition 2.1. Let Xi,... ,Xn be d-dimensional independent random vectors with copula C, such that for b G (0,1) 
known and m = [nb\, the random vectors ACi, ..., Xm have m.c.d.f.s Fi,... ,Fd and the random vectors X^+i, • ■ ■, Xn 
have m.c.d.f.s F[,... ,F'^. 

Then, under Condition \2.1[ the process D„ converges weakly in £“([0, to a stochastic process Dc defined by 

Dc(s,m) = Cys,M) - sCyi,M), (s,m) G [0,l]‘^+\ (13) 

where for (s, u) G [0, 

d 

Cc{s, u) = Zc{s,u) - y] Cj{u)Zc{s, mI-^1) (14) 

1=1 

with Zc is a tight centred Gaussian process with covariance function 

cov{Zc(s, u),Zc(t, n)} = min(s, t){C{u A v) — C(u)C(v)}, (s, u), (t, v) G [0, 

u Av = (min(ui, ui),..., min(itd, Vd)) and itlH = (1,..., 1, ttj, 1,. .., 1). 
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To resample V)n,m, we note that for b < s and u £ [0, l]'^, rewrites as 

'^) “t" '^)} '^n(0; I5 "^) 

and for 5 > s and u £ [0,1]^ as 

H^n,m('55 'W') — Ajt,{s, 1)^71^771,(0, S, 1/-) Aj 7(0, s)'[fC7i^777(s, 6, u) -h Ij '^)}- 

Let B a large integer and consider for (s, t, u) € A x [0,1]^^ and for /3 = 1,..., B, the processes 

[ntj 




nsj +1: [ntJ 


<u)-Ci 


[nsJ +1: \_nt\ 


i= [nsj +1 


,m ('^) ^ 


and 


^n,L{s, t, U) = t,u)- Cj,yns\+l:yut\,rn{u)n^£}^{s, t, 


(15) 


i=i 


with for /3 = 1,..., B, and 1 < i < n, are i.i.d. standard normal random variables. 

Re-sampling versions of ID>„,m can be obtained for /3 = 1,..., B and (s, u) £ [6,1] x [0, l]'^ by 

u) = A„(s, l){Cl%{0, b, u) + ci%{b, s, m)} - A7i(0, s)Cl%(s, 1, u) 
and for (s,u) £ [0,6] x [0,1]^ by 

d 1« (s, u) = A„(s, l)Ci%(0, s, u) - A„(0, s){Cl« (s, 6, u) + Ci%{b, 1,«)}. 

For j = 1,..., d, the functions C'j,[nsj-i-i:[ntj,m appearing in dT51) are an adaptation of the estimator of Cj proposed in 
section 4.2 in lBiicher et al.l (2014) consisting in simple differencing at a bandwidth hk-.i = min{(^ — k + 1)“^^^, 1/2} of the 
empirical copula process: 

Cj,kd,m{u) = -T- Z -, M £ [0,1] , I < k <l <n, 


Y - 


wi th ut = min(u 7 +hk.i, 1) ' Uj = max{uj — hk,i, 0) and = (uj^,... 
of iKoiadinovic et al.l (1201111 . 


,uj, Ud)- This estimator is in spirit of section 3 


Note that is slightly different to 

= A„(s,1)c2'1(0,s,m) - A40,s)c2'1(s,1,m), (s,w) £ [0, l]''+\ 

These resample D^^Jn can be studied in a future research. 

We have the following Proposition (proved in Appendix A.) 

Proposition 2.2. Under the same condition as Provosition \2. 1[ we have the following result: 


(©71,7„, DijL, ■ • ■, -- (Dc, Dg\ ..., D[f)) , 

in £°°(A X [0,1]'^)^’^^, where for {s,t,u) £ A x [0,1]'^, where Dc is defined in (fT^ and D^\...,D^^ are independent 
copies ofCc- 

As a corollary of Proposition 13.II and continuous mapping theorem, we have the following result: 

Corollary 2.1. Let Xi,..., Xn be d-dimensional independent random vectors with copula C, such that for b £ (0,1) 
known and m = [n6J, the random vectors Xi ,..., X^n have c.d.f. F and the random vectors ACtti+i, ..., Xn have c.d.f. 
F'. 


Consider the statistic defined in (HU by 

Sn,m= sup / {]D)77,777 (s,m)}^(7i.„„(m), 

sefo.n Jlo.iU 
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and re-sampling versions of this statistic defined for /3 = 1 ,..., 5 by 


sup 




Under the Condition fOl 

where 


se[o.i] •^[0,1]'^ 

Si%) ^ {Sc, • • ■. 
Sc = sup [ {ID)c(s,m)}^(7(m), 

sGfO.ll 


and Sq \ ..., S^^ are independent copies of Sc- 

An approximate p-value of the test for can be obtained by 


/n”" — 

Pn,B — 


1 f] 1 {si% > S^,m} ■ 


The previous proposition and the Proposition F.l in the supplementary material of lBiicher and Koiadinovid (120151) allow 
to conclude that the test based on Pn,m,B will hold its level asymptotically as n —>■ oo followed by B —>■ oo. 


3. Break detection in the copula when multiple break times in m.c.d.f.s are known 

Suppose in this section that for an integer p > 0 and for j = 1,... + 1, the random vectors Xmj_i+i, ■ • ■, 

have m.c.d.f.s Fij ,..., Fdj where mo = 0, mp+i = n and for j = 1,... ,p mj = lnbj\, 0 < bi <...< bp < 1 known. 
In the sequel, m denotes the vector of break points {nii,... ,mp). Similarly at the section 2, we propose a test for 
Hfp — Ho,c n Hi m where Hi m is defined by 


Hi^m, - j — ,P ^Aij,..., Fdj such that X^j, - - -, ^mj have m.c.d.f.s Fij ,..., F^j . 


Here, for i = 1 ,..., n we consider the random vectors Ufm defined by 

{Flimi ,1 {Xii ), . . . , Fi-mi ^d{^id)') BJq -(-1 — 1 ^ ^ ^ THi 

{Fmi + l:m 2 ,l (-^21)5 ■ • ■ 5 ^ i ^ 1712 

{Fm p + l:n.l {Xii), . . . , Fm p + l-.n,d {X^d)) rup < i < nrip+i = n, 

where for j = 1,... ,p + 1 and q = 1,..., d, Fm^_i+i-.mi,q is the empirical c.d.f. of Xmi_i+iq, ■ • ■, Xm^q defined in ([S]). 
For 1 < k < I < n, denote by Ckd,m the empirical c.d.f. of random vectors Ui^m, ■ ■ ■ ,Un,m- More particularly, for 
1 <k <l <n, l<mi<...< rUp < n and u G [0, l]'^, we have 



Ck-d,m{'^') — 


Ckd{u) 

, 92+1 

IVITtI: 


(m' - m'_i)C'm'_ + 1 ,^' (m) 


mi,...,mp ^ [k,l] 

,..., mq^ e [k, 1] 

m' T —k—1 and m' , t —I 
91-1 92 + 1 


Consider the process defined by 

lll^n,TTi(s, n.) \/nA7j (0, s) (s, 1){C2: [yisj ,771 (^) 1] 


Proposition 3.1. Let Xi,... ,Xn be d-dimensional random vectors with copula C, such that for 0 < bi <...< bp < 1 
known and mj = lnbj\ for j = 1,... ,p, the random vectors Xmj_i+i, • ■ •, Xm^ have m.c.d.f.s Fij ,..., F^j. 


Then, under Condition \2.1l the process Jin.m converges weakly in to Dc defined in (fT3l) . 

The proof of Proposition 13 .1 1 is similar to the the proof of Proposition 12. II in which the supremas are broken on Cp _^_2 
supremas. 
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For s £ [0,1], denote q the integer such that 0 = bo < bi < ... < bg < s < bg+i <...< bp < bp+i = 1. ID)„,rn rewrites 
for (s, u) G [0,1]^^+^ as 


q p 

Dra,m(s, u) = A„ (s,i){Ec n,m {bj-l,bj,u) + Cn,m{bq,S,u)} - A„(0, s){C„,m(s, &g+l, m) + E ^n,m {bj, 6j+i, m)}. 

j=i i=9+i 


For fj = 1,..., B, consider the re-sampling versions 


= A„(s,l){Eci^m(^J-£^i:«) “ An (0, s){C 2 ^^ (s, &,+ i, 16) E , ^J + 1, «)} • 

i=i f=9+i 

Proposition 3.2. Under the same condition as Proposition 3.1, the conclusion of proposition 2.2 holds with'Dn.m instead 
o/]D>n,m and for /3 = 1,... ,B, instead ofOitln. 

This Proposition can be proved in the same way as the prove of Proposition 2.2. 


4. Monte Carlo simulations 


In all the simulations, d-dimensional observations were considered, with either a Clayton (Cl) copula or a Gumbel- 
Hougaard (GH) copula defined for m = (mi, ..., Ud) G [0,1]^ by 


- 1/0 


(n) = max j E Uj — d + 1,0 


0 > 1 , 


Vi=i 


Cf^iu) = exp 


/ 


V 


E{-iog(wj)}® 

f=i 


i/eN 


6» > 0. 


In an equivalent way, the Kendall’s tau of bivariate margins were specified instea d of the parameter 6 of th e copula. 
The Monte Garlo experiments were g enerated using the |R Development Gore TeamI (R statistical system 2013ll and the 
copula package of lHofert et al.l (|20I5^ to sample the Clayton and Gumbel-Hougaard Copulas. The reader may request 
the corresponding routine by contacting the author. 


In a first situation corresponding to the simulations appearing in Table IB. 11 independent samples of sizes n = 
{50,100,200}, and dimensions d = {2,3} are considered, where the first m = [nb\ observations, b = {0.1,0.25,0.5}, 
have for marginal distributions, normal distribntions N(2, 1) and for copula a Clayton copula or a Gumbel-Hougaard cop¬ 
ula. The bivariate margins have a Kendall’s tau of r = {0.25, 0.5, 0.75}. The last n — [n5J observations have for marginal 
distributions, normal distribntions iV(0,1) and for copula a Clayton copula or Gumbel-Hougaard copula and the bivariate 
margins have a Kendall’s tau of r. Typically, the samples were generated under = Hq c H where Bo,c is defined 
in dl]) and Hi^rn in dH]). The percentage of rejection of B™, based on Sn,m defined in (IT^ with a level a = 5% were studied. 


[Table 1 about here.] 

These percentages of rejection are appreciably closed around a = 5%, except for the cases with b = O.I and n = 50,100 
where estimations are calculated on 5 and 10 observations. This may explain the too high percentage of rejection of B™. 

In Table [B^ and IB. 3l independent samples of sizes n = {50,100, 200}, and dimensions d = {2, 3} are considered where 
the first m = [n6J observations, b = {0.1,0.25,0.5} have for marginal distributions, normal distributions fV(2,l) and 
for copula a Clayton copula or Gumbel-Hougaard copula. The bivariate margins have a Kendall’s tau of r = 0.2. The 
last n — \ nb\ observations have for marginal distributions, normal distributions iV(0,1) and for copula Clayton copula or 
Gumbel-Hougaard copula where the bivariate margins have a Kendall’s tau of r = {0.4,0.6}. Typically, the samples were 
generated under alternative hypotheses B™ = Bi^c H Hi^rn where Bi^^ is defined in ([5]) and with 

Hi c ■ There exist k G {1,..., n — 1} and two copulas Ci and C 2 , 

such that Xi ,..., Xk have copula Ci and Xk+i,..., Xn have copula C 2 . (16) 


[Table 2 about here.] 


















[Table 3 about here.] 


The break times k = [ntj, t G {0.1,0.25,0.5} are considered. The percentages of rejection of the hypothesis with 
a level a = 5% are studied. In the same way, the test for Hq — Hq^c H i?o,m where i?o,m is defined in ([3]) is considered, 
based on described in equation ©• The percentages of rejection of based on Sn,m are closed to the percentages of 
rejection of Hq based on S'„. More exactly the percentages of rejection of based on Sn,m are generally smaller than 
percentage of rejection of Hq based for b = (0.1,0.25} and larger for b — 0.5 and t G (0.25, 0.5}. 

Recall that with the hypothesis of a break time m known in the m.c.d.f.s, the rejection of Hq using Sn does not allow 
for a conclusion of a break in the copula of observations contrary to the rejection of using Sn,m- 


5. Discussions and specific situation 

5.1. A Strong mixing condition 


Suppose that the random vectors Xi ,..., are drawn fro m sequences of we akly dependent vectors, in the sense of 
a-mixing dependence (strong mixing dependence) introduced in Rosenblatt ( 195fi[l : 


Definition 5.1. Let (Xi)igz a sequence of random vectors, and for a, 6 G Z = Z U {±c»}, denote by the a-field 
generated by (Xi)a<i<b. The sequence of a-mixing coefficients (orjrsN is defined by 


Or = sup sup \P{Ar\ B) — P{A)P{B)\, r G N. 


The sequence (Xi)igz will be said to be strongly mixing as soon as ar 


r—>-\-cc> 


0 . 


1. The proDOsitions l2.1| and 13.1 [ remain true if we suppose that the marginal probability integral transforms Ui^m, • ■ ■, t/n.m 
defined in are drawn from a strictly stationary sequence {Ui)i^z whose strong mixing coefficients satisfy ar = 
0{n~°‘), a> 1. In this case the covariance structure of Zc is given by cov(Zc(s, m), Zc(t, d)) = min(s, t) cov(l(i7o 
u),l{Uk < v)), (s,m), {t,v) G [0,1]^^+^ 

• 1 Ur),.m. 


2. The propositions l2.2l and 3.2 remains true if we suppose that the marginal probability integral transforms 

defined in ^ are drawn from a strictly stationary sequence {Ui)i^z whose strong mixing coefficients satisfy ar = 
0(n~^ ), a > 3 + 3d/2 and we consider dependent multipliers satisfy (M1)-(M3) appearing in lBiicher and Koiadinovic 
( 20151 section 2) with in = for some 0 < 7 < 1/2 instead of independent multipliers. 


These situations have been stud i ed in Tables IB. 41 and IB. 51 sequences of multipliers were simulated using the pro¬ 
cedure of (|Bucher and Koiadinovid. l2015l The moving average approach. Section 6.1). A standard normal sequence 
of i.i.d. random variables was used in the construction of multipliers. The value of the bandwidth ap pearing in the 


condition (M2) was automatically selected by the proce dure describe d in [[Bucher and Koiadinovid . I2015L Section 5) by 
using the R function bOptEmpProc of npcp package ( Koiadinovid ( 2014ll~l. The ’’comb i ning” function if appearing 


in this same procedure was arbitrarily chosen as if = maximum (see Politis and White . 20041 Section 4). Finally 
the function Lp appearing in the condition (M3) was the convolution product ‘p{x) = Kp * Kp(2x)/Kp * Kp{0), where 
Kp = (1 — 6x^ 6 |a;p)l(|a;| < 1/2) -|- 2(1 — |a;y 1(1/2 < \x\ < 1), a; G K. 


In Table IBAI and Table IBAI dependent samples of sizes n = {100,200} and d = {2,3} are considered with a break 
in the variance of marginal distributions at time m = \nb\, b = {0.1,0.25,0.5}. The samples were generated under 
= Hq^c n i?i,m in Table IB^ and under (-■ilo.c) H in Table [B31 with a break in the copula at time k = [nt\, 
t = {0.1,0.25,0.5}. The data are generated from two autoregressive models (ARl) defined by: 

Xipij = O.bXij £i+ij, j = 1,... ,d and i G Z with Sij white noise. (ARl) 

For i = 1,.. .m and j = 1,... ,d the chosen white noises are N(0, 1), and for z = m + 1,...n the chosen white noises 

are A*(0,16). For b < t, the sample is obtained in the following way: let Ui,... ,Uk +200 be a d-variate i.i.d. sample 

from the copula Ci and i7/c+20i) • ■ • > t/n-i -200 be a d-variate i.i.d. sample from the copula C 2 (under Hff, Ci = C 2 ). 
For i = l,...,m-\- 100, let = ($“^(17a), • ■ •,‘l’~^(C^id)) where $ is the c.d.f. of the standard normal distribution, 
and for z = m + 101,... ,n + 200, Si = (4 x ... ,4 x ^~^{Uid)). Then, Xi = £ 1 , Xm+wi = Em+ioi and for 

z = 1,..., m + 99, and z = m -I- 101,..., n -I- 199 compute recursively 


Xipi — 0.5Aii -|- £ipi. 

Finally, we remove the observations Xi to Xioo and Xm +101 to Xm+ 200 - For b > t, the sample is obtained in the similar 
way. 






































[Table 4 about here.] 
[Table 5 about here.] 


From n = 200, it can seen in Table IB^ that the percentages of rejection of are appreciably closed around a = 5% 
whereas for the whole of break scenarios in copula (Table [B3|, the percentages of rejection of are relatively high. 


5.2. Specific situation 

As an illustration, the bivariate log-returns computed from closing daily quotes of the Dow Jones Industrial Average 
and the Nasdaq Composite for the years 1987 and 1988 have been studied. This is an interesting situation because the 
data highlight a change in the m.c.d.f.s at tim e m = 202 (1987-10-19, corresponding to the ’’Black Monday”). A Cramer- 


von Mises test (see for example iHolmes et al.l . 120131) can allow to confirm this change. Using the procedure described in 
Section 2, an approximate p-value of 0,201 was obtained and no evidence against Hq c is reported. 

Because a marginal gradual change or a multiple marginal change could lead to a rejection of in the case of 
rejection of 77™, the hypothesis of a unique change in marginal distribution should be confirmed. 


5.3. Case of unknown marginal break 

It seem interesting not to fix the break time m = \ nb\ and aggregating this term in the best possible way. Nevertheless, 
such aggregation would be expensive in simulation time, thus the performance of the algorithm should be improved before. 

Another interesting way of a future research will be to consider an estimation of the unknown break time instead of 
m in Sn,m and study the associated statistic. 
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Appendix A. Proof of Proposition [2.11 and Proposition [2T^ 


Proof of Proposition 12.li In the sequel, for fc > ^ the empirical copulas Ck-.i and Cka,m are considered as null by 
convention. Let b G (0,1) such that m = [nb\. 


Le t us consider the two-sided sequential empirical copula process defined in lBiicher and Koiadinovid (|2015ll : lBucher et al.l 
(I 2 OI 4 II by: 

Cnis,t,u) = VnXnis,t){Clns]+l:[nti{u) - C (u)}, is,t,u) S A X [0, 1]'', (A.l) 

where A = {{s,t) G [0, l]^|s < t} . If Xi,..., Xn, are drawn from a i.i.d. sequence (Ai)igz with continuous margins 
Pi ,... ,Fd then (Proposition 3.3 of iBiicher et al. ( 2014 11. 


Pr 

sup |C„(s,<, tt) - C„(s,t, m)|0, 

(s,i,tx)G A X 


with 


C„(s, t, u) = B„(s, t, m) - ^ C'j(M)B„(s, t, 


{s,t,u) e A X [0, ly 


(A.2) 


(A.3) 


and 


ys,t,u) = ^ V {l{U^<u)-C{u)}, 


[nt] 


(s,t,M) e A X [0,1]'^, (A.4) 

where the vector Ui, i = 1 ,..., n, are the vectors of the probability integral tranforms Ui = (Fi(Aii),..., Fd{Xid)). 


Here, we only suppose that Xi, ..., have same copula C. 

Define the process (s,t, it) !->■ Cn,m{s,t,u) similarly to {s,t,u) !->• Cn{s,t,u) in (IA.3I1 and (IA.4() using the vectors 
^[nsj+ 1,7715 • ■ ■ 5 ^[nt\,m instead of ZJ^ns\+i^ ■ ■ ■ 5 ^[nt\ Ir (|A.4p. 

The Proposition 12.11 can be seen as a corollary of the following Lemma: 

Lemma Al. Let Ai,...,A„ be d-dimensional independent random vectors with eopula C, such that for a fixed in¬ 
teger m = \nb\, b G (0,1) known, the random vectors Xi,..., Xm have m.c.d.f.s Fi,..., Fd and the random vectors 
Xm+i, ■ ■ ■, Xn have m.c.d.f.s F[,... ,F'j^. Under Cond,ition \2. 1[ 

sup u') ’Cn,m(s, t/-)| y 0. 

{s,t,u)GA X [0,1]*^ 


consequently, 

Cn,m^Cc zn£°°(Ax [ 0 , 1 ]'^). 

Proof of Lemma EH We will demonstrated the 3 following convergences: 


~ Pr 

sup (s, t, Ia) Cn,mI ^ 0; (-^*^) 

{s,t,u)G{An[ 0 ,b]^)x[ 0 ,l]^ 

" Pr 

sup |^n,m ^n,m ('5; I ^ (-^*^) 

(s,t,'u)G(An[ 6 ,l] 2 )x[ 0 ,l]‘^ 

" Pr 

sup |(Pn,777 (^5 t, It) ^n,m (*^5 L I t 0. (A. 7) 

(s,t,it)G[0,b] X X [0,1]“^ 

On (A n [0, by) X [0, l]'^ or (A 0 [6,1]^) x [0, l]'^, we have C„,m(s, t, u) = C„(s, t, u) and Cn.mis, t, u) = C„(s, t, u), hence 
the desirated convergence follow from (jA.2|) : 


sup |Clyi,77i(s, t, "Ix) (C^ (S, t, tx) I 

(s.t,u)e(An[o.f)]2)x[o.i]‘^ 

= sup |C„(s,t,xi) - C„(s,t, tx)| 

(s,t,ti.)e(An[0,6]2) x [0,1]'^ 

= sup |C;(s,t,M) - C;(s,t,M)| 

(s,t,ii)e{An[0,6]2)x 

< sup |c);(s,t,M) - (C);(s,t,M)| 0 , 

{s,t,u)£Ax[ 0 ,l]^ 


(A.8) 






























where C* and C* are the processes C„ and C„ constructed from Xi ,..., X^nb\ i ^\nb\+i' ■ • ■ > with ..., X* 

such as Xi,..., X* are i.i.d. 

Because by construction, for all {s,t,u) G [0,&] x [6,1] x [0, we have Cn,m(s,t,u) = C„(s, 6 , it) + C„( 6 , t, it) and 
Cn,m(s, t, It) = C„(s, 6 , It) + C„( 6 , t, It), the suprema in (|A.7p is bounded by 

sup |C„(s, 6 ,ii) - C„(s, 6 ,ii)| + sup |C„( 6 ,t,it) - C„( 6 ,t, it)| 

(s,it)e[0,6] X [0,1]*^ (i,it)G[b,l] X [0,1]'^ 

using the triangle inequality, hence the same argumentation as (IA. 8 I) allows us to conclude. 


Using the fact that for any m G {1, ... ,n — 1}, the random vectors Ui^m, ■ • ■, CAi.m are i.i.d., we obtain the weak 
limit of (s.t, It) i- 7 > B„. m.fsA.it) in £°°(A x [0, l]*^) to (s,t,it) !->• Zc(t,it) — Zc(s, it) (see for example the Theorem 
2 .12.1 of van der Vaart and Wellneij . I 2 OOC 1 I I hence the process (s,it) 1 — >■ «) converges weakly in £°°([0, l]'^’*'^) to 

C^(s,it) - sC^(l,it). 


Proof of Proposition 12.21 Proceeding as the proof of Proposition 4.3 in lBiicher et al.l ( 20141 1. for /3 G {1,..., B} and 
(s, t. It) G A X [0, l]'^, put 


lnt\ 

'^u,L{s,t,u) = ^ < m) - C'(lt)}, 


1=1 

and the similar versions and based on [/i,..., [ /„, instead of [/i , Un m- F rom Lemma the fact that 

Ui^m, ■ ■ ■ ,Un,m are i.i.d. and using the Theorem 2.1 in iBiicher and Koiadinovid (1201511 and the continuous mapping 
theorem, we have that 

(c„,™, Ci%, ..., ^ (Cc, cg\ ..., c[f 

in £°°(A X [0,1]'^)^+^. With this result, it is sufficient to demonstrate the following convergences: for /3 = 1,..., i?, 

sup 0n,L{s,b,u)-Cif^l^{s,b,u)\^O, 

(s,ti)G[0,b] X [0,1]^^ 

sup |c£^^( 6 ,s,it)-c 2 ^^( 6 ,s,it)| ^^ 0 , 

(s,ia)G[6,1] X [0,1]^ 

sup |C 2 ^^( 0 ,s,it)-C 2 ^^( 0 ,s,ii)| 0 , sup |C 2 ^^(s, 1 , 11 )-Ci^^(s,l,it)| 0 , 

(s,ia)g[0,&] X [O,!]'^ X [0,1]^^ 

sup |c 2 ^i( 0 , 6 ,it)-c 2 ^^( 0 , 6 ,it)| ^ 0 , sup |(Ci'^^( 6 , l,it)-c 2 ^],,( 6 ,l,it)| ^ 0 . 


In fact, using the same argumentation as (IA. 8 |) and term (B.2) appearing in lBiicher et al.l (I20141) . the previous conver¬ 
gences are automatically verified. 


Appendix B. Simulation study 






















Appendix: simulation study 


Table B.l: Percentage of rejection of based on the Sn,m statistic, computed from 1000 random samples of size n = {50,100, 200} generated 
under H'^, where C is either a d-dimensional Clayton (Cl) or a d-dimensional Gumbel-Hougaard (GH) copula with Kendall’s tau of r. The [nb\ 
first observations have for marginal distributions, normal distributions A^(2, 1) and the n— \nb] last observations have for marginal distributions, 
normal distributions A^(0,1). The test is based on S'^ and replications are computed using independent multipliers 


n 

d 

T 

5 = 0.1 

Cl 

5 = 0.25 

5 = 0.5 

5 = 0.1 

GH 

5 = 0.25 

5 = 0.5 

50 

2 

0.25 

11.3 

9.5 

7.4 

7.9 

6.9 

7.0 



0.50 

16.3 

10.4 

4.8 

10.4 

7.5 

4.0 



0.75 

36.1 

11.3 

5.4 

18.2 

5.8 

1.7 


3 

0.25 

9.3 

9.3 

8.1 

5.3 

6.3 

5.1 



0.50 

10.8 

9.3 

7.0 

3.9 

3.0 

2.5 



0.75 

3.4 

2.1 

1.7 

1.0 

1.3 

0.6 

100 

2 

0.25 

8.2 

7.6 

6.6 

4.3 

2.9 

3.8 



0.50 

8.0 

7.0 

4.9 

6.6 

5.0 

4.4 



0.75 

11.5 

6.0 

2.1 

5.8 

4.1 

1.5 


3 

0.25 

5.7 

6.0 

6.3 

5.4 

5.5 

5.1 



0.50 

8.1 

8.5 

7.9 

3.6 

3.2 

3.4 



0.75 

3.8 

2.7 

1.1 

1.0 

1.2 

0.3 

200 

2 

0.25 

4.7 

5.2 

4.5 

4.9 

5.1 

5.2 



0.50 

4.3 

5.5 

4.2 

5.3 

5.1 

4.4 



0.75 

5.1 

4.0 

1.8 

4.5 

3.1 

1.1 


3 

0.25 

5.2 

5.6 

4.8 

4.0 

4.1 

4.2 



0.50 

7.0 

6.7 

5.5 

3.7 

3.7 

3.1 



0.75 

3.6 

3.2 

2.8 

2.2 

2.0 

1.1 


Table B.2: Percentage of rejection of computed from 1000 random samples of size n = {50,100, 200} generated under Ha = H\^c H 
defined in cs and m, where the first [ntj observations, t G {0.1, 0.25, 05} have for copula a d-dimensional Clayton (Cl) copula with Kendall’s 
tau of 0.2 and the last n — [nt\ have for copula a d-dimensional Clayton copula (Cl) with Kendall’s tau of r. The first m = [nb\ observations, 
b £ {0.1, 0.25,0.5} have for marginal distributions, normal distributions A^(0,1) and the n — m last observations have for marginal distributions, 
normal distributions A^(2,l). Two different tests (based on Sn and S'^) are compared and replications are computed using independent 
multipliers 






Sn,m 

Su 

Sn,m 

Sn 

Sn,m 

Sn 

n 

d 

t 

r 

b = 0.1 

b = 0.1 

6 = 0.25 

b = 0.25 

6 = 0.5 

6 = 0.5 

50 

2 

0.10 

0.4 

15.9 

10.5 

10.5 

15.1 

7.7 

6.7 




0.6 

39.9 

28.5 

23.0 

45.0 

10.7 

8.4 



0.25 

0.4 

21.6 

19.3 

18.0 

21.1 

11.6 

9.6 




0.6 

56.6 

48.5 

48.5 

60.0 

31.8 

19.9 



0.50 

0.4 

23.8 

23.8 

22.0 

25.9 

17.0 

15.3 




0.6 

59.0 

56.3 

55.7 

63.7 

46.3 

43.2 


3 

0.10 

0.4 

12.8 

12.6 

11.3 

15.9 

8.5 

6.0 




0.6 

26.0 

25.4 

20.7 

40.9 

13.6 

8.7 



0.25 

0.4 

20.0 

21.0 

21.7 

25.0 

14.8 

10.5 




0.6 

57.8 

62.6 

56.2 

67.6 

47.2 

30.3 



0.50 

0.4 

28.7 

31.2 

27.5 

35.6 

24.9 

24.6 




0.6 

75.9 

79.4 

76.9 

82.3 

73.5 

72.3 

100 

2 

0.10 

0.4 

10.5 

17.5 

9.5 

27.4 

6.4 

8.1 




0.6 

35.4 

47.9 

28.4 

72.8 

19.5 

25.4 



0.25 

0.4 

23.8 

33.1 

24.0 

40.6 

17.3 

13.6 




0.6 

78.0 

84.2 

74.7 

90.1 

62.7 

47.8 



0.50 

0.4 

32.4 

38.3 

32.4 

45.5 

29.4 

30.0 




0.6 

84.0 

89.7 

85.1 

95.2 

81.9 

85.0 


3 

0.10 

0.4 

11.8 

16.3 

11.5 

22.9 

9.7 

7.9 




0.6 

36.6 

49.4 

33.1 

69.0 

25.7 

21.8 



0.25 

0.4 

30.8 

40.0 

31.1 

45.1 

26.9 

19.7 




0.6 

86.7 

91.9 

87.0 

94.4 

82.4 

63.6 



0.50 

0.4 

43.1 

53.0 

42.4 

59.1 

41.9 

41.8 




0.6 

95.8 

97.0 

95.4 

98.5 

95.6 

96.1 

200 

2 

0.10 

0.4 

11.2 

30.2 

11.7 

45.8 

10.2 

16.6 




0.6 

49.0 

81.4 

42.5 

96.1 

37.7 

70.1 



0.25 

0.4 

35.6 

56.3 

36.2 

68.7 

31.1 

31.7 




0.6 

93.4 

98.0 

93.5 

99.7 

91.1 

89.1 



0.50 

0.4 

50.3 

65.6 

49.8 

76.0 

49.3 

58.4 




0.6 

99.1 

99.9 

98.9 

100.0 

99.0 

99.5 


3 

0.10 

0.4 

12.7 

29.2 

12.6 

44.1 

11.2 

14.3 




0.6 

61.5 

81.6 

59.7 

94.4 

57.9 

59.8 



0.25 

0.4 

46.2 

66.5 

48.8 

72.7 

45.0 

36.7 




0.6 

98.7 

99.9 

99.2 

99.9 

99.0 

94.6 



0.50 

0.4 

69.0 

82.3 

69.3 

88.8 

70.2 

75.2 




0.6 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 











Table B.3: Percentage of rejection of computed from 1000 random samples of size n = {50,100, 200} generated under Ha = H\^c H 
defined in m and m, where the first [ntj observations, t £ {0.1,0.25,05} have for copula a d-dimensional Gumbel-Hougaard (GH) copula 
with Kendall’s tau of 0.2 and the last n — [nt\ have for copula a d-dimensional Gumbel-Hougaard (GH) with Kendall’s tau of r. The first 
m = \nb\ observations, b € {0.1, 0.25, 0.5} have for marginal distributions, normal distributions A^(0,1) and the n — m last observations have 
for marginal distributions, normal distributions A^(2,1). Two different tests (based on Sn and S]^) are compared and replications are computed 
using independent multipliers 






Sn,m 

Sn 

Sn,m 

Sn 

Sn,m 

Sn 

n 

d 

t 

r 

b = 0.1 

b = 0.1 

5 = 0.25 

b = 0.25 

6 = 0.5 

6 = 0.5 

50 

2 

0.10 

0.4 

11.2 

9.0 

8.5 

13.7 

5.6 

4.8 




0.6 

23.5 

22.6 

15.1 

31.4 

8.9 

6.5 



0.25 

0.4 

18.2 

14.9 

15.2 

19.5 

11.1 

7.3 




0.6 

45.9 

41.8 

39.4 

49.1 

28.8 

15.8 



0.50 

0.4 

20.7 

17.9 

19.9 

22.0 

17.1 

15.7 




0.6 

52.9 

55.0 

49.5 

60.1 

46.5 

43.3 


3 

0.10 

0.4 

7.1 

8.2 

6.8 

12.2 

6.1 

5.8 




0.6 

13.4 

16.3 

10.3 

26.1 

8.2 

4.7 



0.25 

0.4 

12.2 

15.8 

14.0 

19.0 

11.5 

8.5 




0.6 

41.3 

49.6 

41.6 

54.5 

36.1 

22.9 



0.50 

0.4 

20.9 

25.9 

21.5 

27.9 

21.8 

18.9 




0.6 

66.1 

72.1 

68.0 

73.8 

63.9 

62.4 

100 

2 

0.10 

0.4 

9.1 

11.5 

7.9 

21.6 

6.6 

6.6 




0.6 

26.5 

44.0 

18.3 

59.5 

14.1 

17.2 



0.25 

0.4 

19.4 

28.3 

20.5 

33.1 

17.1 

11.6 




0.6 

66.7 

76.0 

62.4 

84.2 

55.9 

39.3 



0.50 

0.4 

26.7 

34.6 

27.2 

41.5 

24.8 

25.8 




0.6 

81.6 

87.5 

80.2 

91.0 

79.3 

81.2 


3 

0.10 

0.4 

6.9 

11.9 

6.6 

18.0 

6.2 

5.3 




0.6 

25.1 

38.8 

21.8 

55.1 

19.5 

12.6 



0.25 

0.4 

21.5 

29.4 

22.8 

34.0 

21.2 

14.5 




0.6 

79.1 

87.6 

79.6 

90.0 

76.8 

55.1 



0.50 

0.4 

39.3 

50.8 

38.5 

55.1 

39.3 

39.7 




0.6 

93.2 

95.9 

93.4 

96.9 

92.8 

91.9 

200 

2 

0.10 

0.4 

8.7 

26.4 

8.2 

42.4 

7.8 

15.7 




0.6 

40.5 

79.3 

35.2 

92.7 

33.4 

56.2 



0.25 

0.4 

31.3 

52.1 

29.2 

62.6 

28.3 

26.2 




0.6 

91.1 

98.3 

91.5 

99.8 

90.3 

80.4 



0.50 

0.4 

46.3 

63.0 

45.2 

72.6 

44.2 

51.7 




0.6 

99.1 

99.8 

99.0 

100.0 

99.4 

99.6 


3 

0.10 

0.4 

8.2 

25.3 

9.9 

34.9 

9.6 

9.5 




0.6 

52.5 

80.4 

48.6 

89.4 

48.6 

40.0 



0.25 

0.4 

40.2 

62.2 

42.0 

64.0 

42.0 

31.0 




0.6 

98.0 

99.4 

98.1 

99.6 

97.9 

91.1 



0.50 

0.4 

65.9 

79.5 

65.5 

85.5 

66.1 

69.7 




0.6 

99.9 

100.0 

100.0 

100.0 

100.0 

100.0 











33 m 


Table B.4: Percentage of rejection of computed from 1000 samples of size n = {100, 200} generated under and from two (ARl) models, 
with d-dimensional Clayton (Cl) or Gumbel-Hougaard (GH) stationary copula with Kendall’s tau of r. The m = [nb\ first observations, 
{0.1,0.25,0.5} have stationary margins ^"(0,1) and the n — [nb\ last observations have stationary margins ^"(0,16). The test is based on 
and replications are computed using dependent multipliers 


n 

d 

T 

5 = 0.1 

Cl 

5 = 0.25 

5 = 0.5 

5 = 0.1 

GH 

5 = 0.25 

5 = 0.5 

100 

2 

0.25 

9.0 

10.6 

9.6 


10.3 

8.3 



0.50 

11.2 

8.2 

8.5 


8.1 

6.2 



0.75 

12.2 

11.2 

5.1 


8.1 

5.0 


3 

0.25 

9.7 

9.9 

10.6 

7.0 

10.3 

00 

bo 



0.50 

9.9 

10.7 

8.1 

5.0 

6.8 

7.2 



0.75 

4.5 

3.4 

3.3 

2.5 

2.8 

1.4 

200 

2 

0.25 

5.9 

8.0 

6.5 


6.5 

5.6 



0.50 

6.1 

7.1 

5.1 


5.4 

4.8 



0.75 

4.3 

4.6 

3.4 


2.9 

0.7 


3 

0.25 

6.2 

4.8 

6.4 

5.0 

6.3 

8.2 



0.50 

5.2 

5.7 

4.5 

4.9 

4.3 

4.6 



0.75 

2.1 

3.6 

1.9 

1.0 

1.1 

0.6 









Table B.5: Percentage of rejection of computed from 1000 samples of size n = {100, 200} generated under ->Hq^c H ^i,m s,nd from two 
JARIII models, where the first \nt\ observations, t G {0.1,0.25,05} have for stationary copula a d-dimensional Clayton (Cl) copula (resp. 
Gumbel—Hougaard Copula) with Kendall’s tau of 0.2 and the last n — [nt\ have for stationary copula a bidimensional Clayton copula (Cl) 
(resp. Gumbel—Hougaard Copula) with Kendall’s tau of r. The [nb\ first observations have m.c.d.f. ^(0,1) and the n — [nb\ last observations 
have m.c.d.f. ^"(0,16). The test is based on 5^ and replications are computed using dependent multipliers 

Cl GH 


n 

d 

t 

r 

o 

II 

6 = 0.25 

6 = 0.5 

b = 0.1 

b = 0.25 

b = 0.5 

100 

2 

0.10 

0.4 

14.8 

13.4 

10.5 

12.7 

12.3 

9.1 




0.6 

36.5 

23.9 

10.7 

33.1 

19.6 

9.4 



0.25 

0.4 

24.6 

25.4 

17.5 

22.8 

20.4 

15.4 




0.6 

62.5 

62.8 

45.9 

62.0 

55.7 

45.1 



0.50 

0.4 

31.9 

32.2 

25.2 

27.0 

24.8 

25.8 




0.6 

74.7 

71.0 

69.1 

73.8 

72.1 

66.4 


3 

0.10 

0.4 

12.2 

13.2 

11.8 

11.3 

10.6 

9.0 




0.6 

34.9 

23.1 

13.8 

26.0 

18.3 

11.1 



0.25 

0.4 

23.7 

31.5 

20.9 

23.7 

25.3 

18.7 




0.6 

74.0 

77.0 

64.9 

67.3 

72.6 

59.0 



0.50 

0.4 

38.8 

40.8 

36.1 

32.3 

35.0 

35.5 




0.6 

85.1 

87.6 

86.8 

85.6 

86.6 

84.8 

200 

2 

0.10 

0.4 

10.8 

10.4 

9.6 

10.5 

8.6 

7.9 




0.6 

39.9 

21.9 

11.3 

36.3 

20.1 

12.2 



0.25 

0.4 

27.9 

29.5 

22.1 

24.2 

26.6 

20.6 




0.6 

80.9 

81.5 

72.1 

79.8 

82.0 

71.9 



0.50 

0.4 

39.5 

37.5 

35.7 

35.7 

38.4 

32.6 




0.6 

91.9 

92.5 

90.3 

92.5 

90.9 

92.0 


3 

0.10 

0.4 

14.9 

11.4 

8.9 

10.4 

9.6 

8.4 




0.6 

45.8 

26.6 

17.6 

45.2 

21.4 

15.6 



0.25 

0.4 

37.3 

38.9 

31.9 

30.6 

35.0 

29.1 




0.6 

93.5 

94.1 

86.8 

90.7 

93.6 

86.9 



0.50 

0.4 

54.6 

53.0 

54.0 

51.5 

51.8 

51.9 




0.6 

99.3 

98.0 

98.0 

98.5 

98.9 

98.4 













